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Abstract 

We present compact, analytic expressions for the age-redshift relation t(z) for standard 
Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology. The new expressions are given in 
terms of incomplete Legendre elliptic integrals and evaluate much faster than by direct numer- 
ical integration. 
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I. INTRODUCTION 



Since Type la supernova observations [Ml [19| have favored a Universe with A / 0, 



interest in FLRW cosmologies has grown. For this family of models most observational 
relations, e.g., the Hubble curve, are given by integral expressions; however, Kantowski et 
al. || recently succeeded in giving useful analytic distance-redshift relations for them. In 
this paper we extend these analytic results to include the age-redshift relation t(z). These 
new expressions are useful for any computation that requires a transformation t(z) from 
the observed redshift variable z to the age r of the Universe at that z. Feige |7[] provides 
related light travel times in terms of Legendre elliptic integrals; however, his expressions are 
not easy to make use of. A presentation closer to what we give appears in ||. There, light 
travel time as a function of redshift was given for the Qq = 1 universe, see ([18]). For all other 
cases, Edwards § was only able to give r and z parametrically as Jacobi elliptic functions 
of conformal time u = J dt/R{t). 1 

We have concentrated on giving useful and valid expressions for r(z) appropriate for 
all big bang models in the first quadrant of the Q m -Q\ plane. Because the incomplete 
Legendre elliptic integrals have branch points, more than one expression is necessary to 
completely cover this domain, e.g., see (H) and fllPP. These new expressions should be quite 



useful for everything from gravitational lensing to high z evolution studies. As an example 
|T3| , |TT| , [L8| and || all estimate event rates of supernovae at very high (z > 1) redshifts, 
given an observed star formation rate. Computing such event rates for any choice of (Q m ,flA) 
requires the transformation r(z). We have made similar estimates of event rates and find 
that our computations are reduced from hours down to minutes when our analytic t(z) is 
used. 2 In §2 we present our results and in §3 some conclusions. 



1 The error detected when numerical checks were performed by Campusano et al. Q was evidently caused 
by an error in equation 361.54 of Byrd & Friedman |2j, see footnote 3 of Q 

2 FORTRAN 90 and Mathematica implementations of the results presented here are available at 



http: / / www. nhn.ou.edu/^thomas/z2t. html 
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II. AGE OF THE UNIVERSE IN TERMS OF LEGENDRE ELLIPTIC INTE- 
GRALS 



The expression for the age of the Universe at the time a source at redshift z emits light 



is 



1 f°° 

T(Q m ,flA;z) = — I — 



dz 



(1 + 2)^(1 + zf{l + n m z) - z(2 + z)Q A ' 
and can easily be derived. For a Friedmann- Robertson- Walker (FRW) universe, i.e., Q\ 
(0) can be integrated in terms of elementary functions, 
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This expression is also valid for the Einstein-de Sitter universe, i.e., limit Q m — > 1, as well 
as when Q m > 1. For massless big bang models, fl m = and < < 1, the integral is, 
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When Q\ ^ and Q m ^ 0, becomes an incomplete elliptic integral and can at best be ex- 
pressed as a combination of the three independent Legendre elliptic integrals F((p, k), E(cf), k), 
and n(0, a 2 , k). 3 The form of the resulting expression depends on what portion of the fl m - 
fl\ plane is being investigated. Because the cubic under the radical in ([!]) is the same as that 



contained in integrals for the luminosity distance as given by || and |10| a similar analysis 
is required. Below we outline results, hoping to make our expressions easy to use. 

a. n = n m + n A / 1 

As seen below, t(z) depends on (fl m , Q\) primarily through a single parameter b 



-(27/2). 



n 2 n A 

in, A 



(4) 



;i-o ) 3 

This parameter divides the Q m -Q\ plane (see Fig. 1) into four domains where the results of 
integrating ([!]) differ. We will ignore one of the four domains and its b = 2 boundary where 



Only two of the three are needed to give t(z) and they are defined by: F((f>, k) = J Q 1/yl-k 2 sin cf> d(f> 
and 11(0, a 2 , k) = L 1/ (1 — a 2 sin 2 (f>) \f\ — k 2 sin 2 4> d<fi. The particular integrals needed can be found 
in §. 
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big bangs don't occur. In the following we use the familiar parameter k = (Qq — l)/|fio — 1 1 , 
which is determined by the sign of the 3-curvature, to distinguish between open and closed 
models. When b < 0, k = — 1 and when b > 0, k = +1. Results for special boundaries other 
than f^A = 0, i.e., (Q) and Q m = 0, i.e., ([$]), are given in subsection B below. The three 
special boundaries needed are: b = ±oo i.e., Q = Q m + Q A = 1; b = 2; and b = 27(2 + y / 2)/8. 

1. Results for the two large domains, b < and 2 < b, can be combined by defining 
intermediate constants v K , y\ and A: 

^-l) + v /fe(6-2)) 1/3 , (5) 
2/i = o > ( 6 ) 



A =v/yi(3 yi + 2). (7) 

These constants depend on 6 alone and are only used for convenience of presentation. If the 
reader desires, the source of these parameters can be found in ||10|| . For this case we give two 
expressions for the integral ([[]). Both are valid except for special combinations of Q m ,Q^, 
and z. If one fails the other can be used. These expressions fail when 11(0, a 2 , k) and the 
logarithm have canceling infinities. Both can simultaneously fail only when Q m , Q\ values 
are on the b = 27(2 + v / 2)/8 curve and then only for a specific value of z (see Figs. 2 and 
3). This special b case is given in ( p2|) and is good for any z. The first expression is: 
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where 
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The second expression is obtained from the first by using a "special addition formula" 
analytically extended from a corrected version of equation 17.03 of 0. This transformation 
changes the a 2 value of 11(0, a 2 , k) and hence moves the associated branch point. The 
resulting second expression is: 
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respectively the modulus and argument of the elliptic integrals, are 
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In (|8p and (|TUp the z dependence of r is contained in <p z , hj, and /i z . All other terms 
depend on Q m and ft\, and are easily evaluated using (|)-(0). In Figures 2 and 3 the dotted 
lines show points where the first expression @ fails for z = 1 and z = 2 respectively. Failure 
of the second expression ( |lOl) is shown by the dashed lines. Notice that these curves always 
intersect somewhere on the b = 27(2 + y/2)/8 curve for a common redshift. 

2. If < b < 2, we define the three different intermediate parameters j/j, ?/2 an d 2/3 
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For this case ([I]) integrates to give 



where k and 



are defined by 
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V yi-(l + z)O m /(l-lV)' (17) 
In ( |15|) the 2; dependence of r is contained in <f> z . All other terms depend on Q m and fi\., 
and are easily evaluated using (|) and ([[4]). There are two domains in the Q m -Q^ plane 
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where < b < 2; however, the result for this case fllq) applies only to those models which 
have big bangs. 
B. Special Cases 

1. fig = f2 m + fiy\ = 1 

This is the spatially flat model (b — > ±oo) and for it the age-redshift integral takes on 
a simpler form. This result is easily obtained using elementary integration methods. This 
result is well known: 
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2. b = 2 

This value of 6 can be identified with "critical" values of the cosmic parameters [H]. The 



following result is equivalent to the 6 = 2 value given in (|15D ; however, it is a much simpler 
expression: 
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This r(z) doesn't apply to the Einstein-Lemaitre universe (b = 2 where > fi m /2) [p~2| , 
which starts expanding from the finite static Einstein radius at t — — oo. However, it does 
apply to the fl\ < fl m /2 models which start with a big bang and expand to the Einstein 
radius at t — +oo. In the Q m -Q\ plane the static Einstein universe itself is a point at oo 
on the two 6 = 2 curves where Q m /Q\ — > 2. If wanted, the b = constant > 2 curves can be 
drawn using the following expressions. Because fl\(fl m ) is double valued, two expressions 
must be given. For the upper part of the curve: 



fi — 1 = 3a/2/& Vt m cosh 



cosh 



6/2 (o- 1 - 1; 



(20) 



where < Vt m < 1/(1 — ^2/b). In this expression hyperbolic cosine analytically becomes 
cosine for Q m > 1/(1 + y2/b). 
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For the lower part of the curve: 

fio — 1 = 3y2/6 tt m cos 
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where 1 < Q m < 1/(1 — y2/b). For 6 = 2 (see Fig. 1) the max value of Q m is 'oo' (the static 
Einstein universe); however, for the next case (see Fig. 3) the upper and lower parts of the 
curve meet at finite Q m ~ 1.7. 
3. 6 = 27(2 + V2)/8 

This result is equivalent to the values given by @ and (|Tt]) except for certain redshifts 
where canceling infinities appear in 11(0, a 2 , k) and the respective logarithms. It is a simpler 
expression and is valid for all z values, for this particular 6, 
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We found it necessary to compute r(z) for this particular 6 value to overcome the occasional 
simultaneous failures of (|j) and flTDD, see Figures 2 and 3. 



III. CONCLUSIONS 

We have given valid analytic expressions for t(z) in FLRW, the age of the Universe as a 
function of redshift, which are relatively simple and are quite useful when a fast computer 
implementation is needed. 4 These expressions completely cover the big bang models of the 
first quadrant of the (fi m , fi \) plane. If lookback times are wanted they can additionally 



FORTRAN 90 implementation of the results presented here and available at 



http://www.nhn.ou.edu/~thomas/z2t.html are 20-40 times faster than a traditional Bulirsch-Stoer 



integrator |15|. 



7 



be obtained from results given here by simply evaluating r(0) — t(z). Readers that are 
interested in adding radiation pressure as a source of gravity should see [0] and || and cited 
references. 

Even though we give several expressions for t(z), most of the (Q m ,Q^) plane, which 
includes currently favored values, is covered by case Al, i.e., result (Q) or flTDp. If r(z) for 
the flat model, Qq = 1, is wanted, the simpler result flUf) should be used. Results for Q ^ 1, 



(I), ([TOD , and (|T^), appear complicated because of the presence of extra constants, e.g., A 
and yi that have been retained to compactify formulas. The reader should keep in mind 
that these are simply constants that depend on (fl m , ^a) through the single combination b 
of ([|). We could have eliminated these auxiliary constants and given t(z) directly in terms 
of the two parameters fl m and Q\] however, such expressions would take up more than a 
page. 

Expressions (|8|) and (|T0|) for t(z) remain real but as presented can contain imaginary 
terms because of branch points. The threshold is defined by 1 — a 2 sin <f) z = Oinn(0 2 ,a 2 ,k). 
If 1 — a 2 sin 2 cj) z < canceling imaginary terms appear in IT(0 2 , a 2 , k) and the logarithm. For 
expressions that avoid this imaginary complication the reader simply replaces the argument 
of the logarithm with its magnitude and U(<f) z ,a 2 ,k) with its principal part. At threshold 
points where 1 — a 2 sin 2 2 = 0, canceling infinities appear in IT(0 2 , a 2 , k) and the logarithm. 
The infinity problem is avoided by switching between (^) and (0). If both have infinities 
then gives the correct result. 
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FIG. 1: The £l m -Q\ plane showing various b domains that require different expressions for age- 



redshift t(z) for standard FLRW. Expressions (|S[) and ( |10|) are both appropriate for b < and 
b > 2, while ( |i"8| ) is appropriate for < b < 2. Simpler expressions exist for various boundaries: 
n A = (0), n m = (§), 6 oo O O = 1 (0), and 6 = 2 (0). 
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m 

FIG. 2: The £l m —Q\ plane showing curves where age-redshift t(z) expressions (g) and ( |i~0| ) fail 
for redshift z = 1. Equation (||) fails along the dotted curve and (|l(]) fails along the dashed curve. 
Both fail where they intersect on the b = 27(2 + y/2)/8 curve; however, ( [22] ) gives the t(z) value 
at any point on this curve for all redshifts. 
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FIG. 3: Same as Fig. 2 except for redshift z = 2. 



